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We investigate the power-dependent photoluminescence spectra from a strongly coupled quantum 
dot-cavity system using a quantum master equation technique that accounts for incoherent pumping, 
pure dephasing, and fermion or boson statistics. Analytical spectra at the one-photon correlation 
level and the numerically exact multi-photon spectra for fermions are presented. We compare to 
recent experiments on a quantum dot-micropiller cavity system and show that an excellent fit to 
the data can be obtained by varying only the incoherent pump rates in direct correspondence with 
the experiments. Our theory and experiments together show a clear and systematic way of studying 
stimulated-emission induced broadening and anharmonic cavity-QED. 
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Introduction.- Single quantum dot (QD) - cavity sys- 
tems facilitate the realization of solid state qubits (quan- 
tum bits) and have applications for producing single 
photons [l|, [2I, Q and entangled photons 0, [5[. Rich 
in physics and potential applications, the coupled QD- 
cavity has been inspiring theoretical and experimental 
groups to probe deeper into the underlying physics of 
both weak and strong coupling regimes of semiconduc- 
tor cavity-QED (quantum electrodynamics). Key sig- 
natures of cavity-QED include the Purcell effect and 
vacuum Rabi oscillations. Although a well known phe- 
nomenon in atomic cavity optics [6[ , vacuum Rabi split- 
ting in a semiconductor cavity was only realized a few 
years ago 0, B 0]- Inspired by the recent surge of re- 
lated experiments, many researchers have been working 
hard to develop new theoretical tools to understand the 
semiconductor cavity-QED systems. For example, the 
persistent excitation of the cavity mode for large exciton- 
cavity detunings was measured [10, [U , and qualitatively 
explained by extended theoretical approaches that ac- 
count for coupling between the leaky cavity mode and the 
exciton, and by showing that the main contribution to 
the emitted spectrum comes from the cavity-mode emis- 
sion [12I, |T^, 0, [15I, llB] . These formalisms assume an ini- 
tially excited exciton or an initially excited leaky cavity 
mode, and they are valid for low pump powers. However, 
an interes ting question that has been posed recently, e.g., 
see Refs. [Ulillii, is what is the role of an incoher- 
ent pump on the photoluminescence (PL) spectra, where 
the pump can excite the exciton or cavity mode? To 
experimentally investigate the pump-dependent spectra, 
two recent experiments have been respectively reported 
by Miinch et al. )2p] for a QD-micropillar system, and 
by Laucht et al. |2l| for a QD-photonic crystal system; 
these measurements show the pump-induced crossover 
from strong to weak coupling. 

In this work, we present a master equation (ME) the- 



ory that self-consistently includes incoherent pumping, 
stimulated emission, and pure dephasing. We derive an- 
alytical results at the level of one-photon correlations and 
present numerically exact results for the multi-photon 
spectra. We analyze the Wiirzburg (2^ experiments di- 
rectly and show the striking differences with previous 
models that neglect the direct influence of stimulated 
emission [l8|, [l9|, |2l|. For the incoherent pumping of 
the exciton, we present two ME models: a thermal bath 
model (c.f. a two-level system) and a heat bath model at 
large negative temperature (c.f. a multi- level laser sys- 
tem). Accounting for fermion statistics, pure dephasing, 
and the thermal bath model, an excellent fit to the data is 
obtained by only changing the incoherent pump rates in 
direct correspondence with the experiments. 



Cavity system and model.- The system investigated 
here is shown as a scanning electron microscope (SEM) 
image in Fig.^ along with the extended experimental 
data of Ref. [20]. We make the following assumptions: 
the cavity is single-mode in the frequency of interest; the 
coupling between the cavity and target QD exciton is de- 
scribed through a coupling rate g; the decay rate of cav- 
ity is Fc; for the strongly coupled QD, we include only 
the target exciton as a system operator, and consider 
both radiative decay, F^^, and pure dephasing, F^. The 
QD-cavity system is driven simultaneously by an exciton 
pump, Px^ and a cavity pump, Pc] the former is caused 
by the incoherent relaxation of electron-hole pairs from 
the higher energy level, and the latter is due to the cav- 
ity coupling with off-resonant excitons (probably coming 
from other QDs in the cavity layer). To treat the inco- 
herent excitation, we consider a system-reservoir interac- 
tion [m , apply a Born-Markov approximation, and trace 
over the cavity and target exciton pump reservoirs {bath 
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FIG. 1: (Color online) Typical broadband PL spectrum that 
is emitted when a target exciton is closely resonant with the 
cavity mode (near coq = 1331.355 meV); away from the target 
exciton, there are a series of other exciton levels that can also 
couple, off- resonantly, to the cavity mode. The SEM image 
shows our micropillar cavity and the QD layer. The emitted 
photons from the QDs are detected through vertical emission. 



Pi2 = Px and P21 =^x-^Px'', model-2 (heat bath at large 
negative temperatures) [24]: P12 = Px and P21 = T^; 
Px is the target exciton pump rate which is presumed to 
proportionally follow the experimental pump power. 

One can next derive analytical spectra at the level 
of one-photon correlations, or compute the exact nu- 
merical s pec tra for n— photon correlations, e.g., see 
Refs. [H, HJ). We will present both approaches. Us- 
ing Eq. (2), adopting the one photon-correlation approx- 
imation {cFzo) = — (a), and applying fermion statis- 
tics = 1, we exploit the quantum regres- 
sion theorem [22[ to derive the equation of motion for 
the two-time correlation functions, d {0) {t)a{t ^ r)) / dr 
and d {d^ {t)a~ {t -\- r)) /dr. Subsequently, the steady- 
state form of the dominant cavity-emitted spectrum [1^ 
is obtained from 5'cav(^,^) = Pcav(^) 'S'cav(^), with 
^eav(c^) = r,/^lim,^ooRe{/o"(at(t)a(t + r))e^--dr}, 
where Pcav(^) is a geometrical factor that depends on 
the detector /collection optics. One obtains 



5'cav(^) = — Re 

TT 



C{uo)D{uo) -g^ C{uo)D{uo) - 



,(3) 



approximation). We have 



dp 
~dt 



n 



(1) 



with the 

hujcd^d + 



system Hamiltonian, Hg 



huOxCF^CF 



hg{a d^ + (j+a), where d represents the cav- 



ity mode operator, a+/ are the Pauli operators of the 
target QD exciton (with resonance frequency cox)^ and 
ujc is the eigenfrequency of the leaky cavity mode. The 
target exciton and cavity mode get pumped incoherently 
through the corresponding reservoirs. The state of the 
reservoirs can be written as pp = J2k Pkk l^fc) (^fcl' 
O = X, c; where is the density of reservoir modes 
and is number of photons in the mode of wave vec- 
tor k. The correlations for the photon reservoir opera- 
tors dl are given by (a^) = 0, {{dD^dl,) = nlSkk' , and 
(a^(a/e)'*') = (n^ + l)^/c/e'- Defining the average pump 
photon number around the cavity frequency as = n^, 
at = cjc/c, yields the effective incoherent cavity pump 
rate: Pc = Fcn^. This incoherent pump process agrees 
with the model of Tian and Carmichael [23] • The super- 
operator in Eq. (1) becomes 



C{p) = {2d^ pd — dd^ p — pdd^) , 

- ^ {2apd^ — d^dp — pd^d) , 

- (2a~^pa~ — a~a~^p — pa~a~^) 

- ^ {2a- pa^ - a^a-p - pa^a') + ^ (a.pa, -p) , (2) 



which is in Lindblad form. For the exciton pump we 
consider two different models: model- 1 (thermal bath): 



where C{u;) = uo — uOc^^Tc^ and D{uo) = oo — uOx^^{P2i + 
P12 +r^). The subscript 'ss' represents the steady-state 
solutions, that are given by 



= \- I (4) 

g^T{P2l+Pl2 +re) +re(P21+Pl2)(x + Acj' 



_ -i5({ata),, - -p£^^) (iA 



cx + 2, 



r2 



4 ' ^CX ^ P21+P12 

Pi2 + i9{{a^a-),,-{aa+)J 
P21 + P12 



(5) 
(6) 



where T = P21 + P12 + + Tc, and Acx = - ^x- To 
recover boson statistics, one simply replaces the P21 +P12 
terms above by P21 — P12 and sets T^ to zero; We 
stress that the above formulas are substantially differ- 
ent to previous models that neglect stimulated emis- 
sion [18, 19, 21]; in particular, we have no unphysical be- 
havior as Tc = Pc, and we get qualitatively different sat- 
uration behavior of the QD exciton. Similar incoherent 
pump models, with pump-induced stimulated emission, 
have also been proposed recently by Ridolfo et al. [2^, 
though they concentrate exclusively on the model-2 ex- 
citon pump and they neglect pure dephasing; thus their 
analytical formula applies only to a boson system. 

The power dependent PL.- To highlight the under- 
lying physics of pump-induced PL, we proportionally 
change Px (and Pc) in our model, and keep all other 
parameters fixed (i.e., ^, T^^ T^). The fixed parame- 
ters are either known for our experimental system, e.g. 
Fx = 0.002 meV [25|, or are accurately obtained from the 
fitting the experimental data at low powers, where g = 
0.045 meV and Tc = 0.08 meV. We have also included a 
dominant pure dephasing exciton decay, T^ = 0.035 meV, 



3 




-0.2 0.2 -0.2 0.2 
(o-coq (meV) 




/ 2 




2 

P/g 



1331.4 1331.55 
CO (meV) 



FIG. 2: (Color online) The on-resonance (ujc ~ ^x) PL 
spectra, for different excitation powers, (a) Solution of our 
ME with model 1 (left) and model 2 (right). The red 
curve is the one photon spectra and the blue curve is the 
multi-photon case. The bottom-to-top panels have Px — 
[0.12, 0.5, 4, 16, 64] 0.02125^ (0.0003 - 1.36^), and Pc = 1.6P^. 
(b) ME solution without stimulated emission, (c) Mean exci- 
ton number (dashed) and photon number (solid), for our ME 
with model 1 (left); ME solution without stimulated emission 
using model 1 (right), (d) Experimental data corresponding 
to Pexp = [0.12,0.25,0.5,2,4,8, 16, 32,64] /iW, and model-1 
fits (multi-photon and stimulated emission included), where 
Px/c proportionally follows the experimental values; the inset 
shows the integrated PL (experiment and theory). 



caused by electron-phonon scattering and spectral diffu- 
sion. The chosen values of Px range from 0.003 — 1.36 
and Pc = l.QPx' The justification for allowing Pc to also 
follow the power of the laser is due to the fact that our 
micropillar measurements show a clear linear dependence 
with power for the cavity mode. For other QD-cavity sys- 
tems, such as for a few QDs in a photonic crystal cavity, 
Pc may saturate at much lower powers. In Fig. 2(a), 
we first show the power-dependent spectra for model-1 
(left) and model-2 (right); and in Fig. 2(b), we compare 
the trend expected from a ME model that neglects stim- 
ulated emission processes. The red curves show the one 
photon results and the blue curves show the multi-photon 
case. Although all figures show a similar trend of the dou- 
blet becoming a singlet as a function of power, the high 
power linewidths are substantially different. In particu- 
lar, the model with stimulated emission predicts a much 
larger pump-induced broadening as a function of power. 
In the absence of stimulated emission, the pump-induced 
broadening is suppressed, and the larger pump rates re- 
sult in negative exciton and photon densities. The mean 
exciton number (dashed) and photon number (solid) are 
shown in Fig. 2(c) using the multi-photon model. Here 
we see the drastic influence on the predicted densities if 
stimulated emission is not included (right), where nega- 
tive photon densities are predicted in addition to regimes 
oiux > 1, both of which are obviously unphysical; though 
we model-1 densities, model-2 gives similar unphysical re- 
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FIG. 3: (color online) The on-resonance (cjc ~ ^x) PL spec- 
tra, for different excitation powers, but for a boson model, 
(a) ME with stimulated emission, using exciton pump model- 
1 (blue) and model-2 (red), (c) Corresponding mean density 
plots: exciton number (dashed) and photon numbers (solid); 
for clarity the model-1 densities are multiplied by 100. (b,d) 
As in (a,c), but without stimulated emission. 



suits [27|. Of course, with stimulated emission neglected 
in the model, the regime of Pc > Tc is phenomenolog- 
ically not allowed [1^, so the top spectra in Fig. 2(b) 
are not reliable. It is interesting to note that, even for 
pump rates as small at Px = 0.085^, multi-photon states 
(c.f. Jaynes-Cummings model [22]) are already important 
contributions to the nonlinear spectra, and we find that 
2-4 photon states are enough to get good convergence. 

The experimental data is shown in Fig. 2(d), along- 
side our fermion model-1, and there is an excellent cor- 
respondence. We stress that the only fitting parameter 
is a proportionality constant. Although may also be 
pump-dependent, we find that increasing its value by 1-2 
orders of magnitude has little influence on our high-power 
PL, as the stimulated-emission-induced broadening is by 
far the dominant source of broadening. To have further 
confidence in the theory, it is important that the models 
consistently fit the normalized PL, on and off-resonance, 
as well as the integrated PL. We obtain very good fits to 
the spectra when the cavity and exciton are off resonance 
(not shown) and the integrated PL [shown as an inset in 
Fig 2(d)], without changing any parameters. 

Since our QDs are rather large, e.g., elongated with 
lengths ^n the order of 100 nm and widths of about 
30 nm [7], it is natural to present the nonlinear boson 
PL calculations as well. In Fig. 3 we display the ex- 
act boson PL using exciton-pump models 1 and 2, again 
with and without stimulated emission terms. Since pure 
dephasing cannot be included, we set Tx F^^ + F^. 
Clearly, none of the PL follow the trends of the exper- 
iments, and only the thermal bath models produce net 
positive densities for all pump rates. Moreover, even the 
low PL have different lineshapes due to the important 
effect of pure dephasing, which acts to suppress the Rabi 
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FIG. 4: (color online) (a) Mean exciton number versus Px 
(with Pc = 1.6Px, as before): exciton pump model- 1 (solid) 
and model-2 (dashed), (b) Corresponding mean photon num- 
ber (left axis: blue) and Fano factor F (right axis: red); 

oscillations without affecting the envelope of the popula- 
tion decay. While it has been discussed before that the 
boson model (with model-2) [20j apparently fits well to 
the same data under variation of the coupling constant g 
and three other free parameters {Tx-, Px-, Pc)] we believe 
that having so many free parameters can be detrimen- 
tal to highlighting the underlying physics. We conclude 
that our nonlinear PL spectra unambiguously follow the 
presented fermion model, and we are thus well into the 
regime of anharmonic cavity QED. 

High pump-power inversion and lasing- Finally, we 
briefly connect to the prospects for observing one exci- 
ton lasing in such a QD system. It Js well known in the 
field of atomic optics, e.g., see Ref. [2^, that the spectral 
properties of pump-dependent PL can be investigated to 
explore the regime of single atom lasing. Characteris- 
tic signatures of single state lasing in atomic physics in- 
clude spectral narrowing, inversion, and a regime of lin- 
early increasing mean photon number as a function of 
pump power. On the other hand, an incoherent pump 
of thermal photons will naturally be detrimental to the 
prospect of achieving single photon lasing. In Fig. 4(a) we 
use model- 1 (solid) and model-2 (dashed) to investigate 
the pump-dependent mean exciton number and the mean 



photon number (panel (b)). As expected model- 1 (ther- 
mal bath) does not exhibit any inversion, though model-2 
does allow inversion. Both models allow a mean photon 
number of greater than 1. However, the Fano function 
(photon number variance) [28] shows no evidence of a 
maximum, and thus there is no lasing threshold in this 
system. To achieve single exciton lasing with the present 
model/system, we have numerically verified that one re- 
quires a much smaller Pc/Px ratio and a significantly 
smaller Fc; for example, Pc = and Fc = 0.01 meV gives 
a very clear lasing threshold and order-of-magnitude re- 
ductions in the PL linewidth. Experimental activity on 
single QD lasers has begun and, in future work, we 
will explore the key signatures of single exciton lasing 
using a more detailed multi-level excitation scheme. 

Conclusions.- A master equation formalism, with in- 
coherent pumping, pure dephasing, and a QD fermion 
model, has been introduced and used to investigate the 
power-dependent PL spectrum of a QD exciton under 
steady-state pumping. We have shown the importance of 
self-consistently including stimulated emission, and vali- 
dated our model by directly comparing with recent exper- 
imental data on semiconductor micropillar-cavities. Us- 
ing the proposed thermal bath model, an excellent fit to 
the data is obtained by only changing the pump rates 
in direct correspondence with the experiments, showing 
that we are well into the elusive regime of anharmonic 
cavity- QED. Moreover, we have shown that our exci- 
tation models produce positive-definite densities for all 
pump rates. 
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